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The interactions of a i5/yl=70-MeV beam with a Be target was used to populate levels in Be, 
B, and C isotopes which undergo decay into many-particle exit channels. The decay products were 
detected in the HiRA array and the level energies were identihed hom their invariant mass. Corre- 
lations between the decay products were examined to deduce the nature of the decays, specifically 
to what extent all the fragments were created in one prompt step or whether the disintegration pro- 
ceeded in a sequential fashion through long-lived intermediate states. In the latter case, information 
on the spin of the level was also obtained. Of particular interest is the 5-body decay of the *C 
ground state which was found to disintegrate in two steps of two-proton decay passing through the 
''Beg.s. intermediate state. The isobaric analog of ^Cg.a. in ®B was also found to undergo two-proton 
decay to the isobaric analog of ^Beg.s. in «Li. A 9.69-MeV state in ^^C was found to undergo prompt 
4-body decay to the 2p~f 2q exit channel. The two protons were found to have a strong enhancement 
in the diproton region and the relative energies of all four p-a pairs were consistent with the ^Lig.s. 
resonance. 

PACS numbers: 21.10.-k,25.70.Ef,25.60.-t,27.20.-|-n 
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I. INTRODUCTION 

The ground and excited states of many light proton- 
rich nuclei decay by the emission of protons and other 
charged particles. A number of these states disintegrate 
into just protons and alpha particles and other very light 
fragments. Well known examples are the ground states 
of ^Be and ^B which decay to the 2p-\-a and the p+2a 
exit channels. Such states can be studied with resonance 
decay spectroscopy where all of the decay fragments are 
detected. The excitation energy of the parent state can 
be determined from the invariant mass of the detected 
decay products. For exit channels with more than two 
fragments, the angular and energy correlations between 
the fragments can provide information on the spin of the 
state and the nature of its decay, i.e., whether the de- 
cay proceeds through a series of sequential 2-body decay 
steps or via a more prompt process. 

In this work we have utilized a ^C beam to study a 
number of such states which were formed via nucleon 
knockout and more complicated reactions. Of particular 
interest is the ground state of *C which is known to be 
unstable to the 5-body exit channel Ap+a. Information 



about ^C is rather sparse and the 5-body exit channel 
has never been observed experimentally. A particular 
question for ^C disintegration and other many-particle 
exit channels detected in this work is whether the final 
particles were created in one fast process or in a series 
of sequential steps with the creation of long-lived inter- 
mediate states? Of course there is not always a clear 
demarcation between such processes. If the lifetime of 
the intermediate states becomes too short, does it make 
sense to still describe the decay as two or more sequential 
steps? 

One situation that is particularly clear is the two- 
proton decay described by Goldansky Tj| where no inter- 
mediate state can be energetically accessed. Examples of 
this are the two-proton decays of ^^Fe [3| and ^^Zn^Si]. In 
cases where this condition is not reached, Bochkarev et 
al. considered a distinction in the decay type based on 
the width of the potential intermediate state. For a 3- 
body exit channel, if the potential intermediate state has 
a width T I which is of the same magnitude as the kinetic 
energy E], released in the decay to this state, then its de- 
cay was deemed democratic, i.e., all energy scales in the 
subsystems and the three-body systems are comparable. 
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We have constructed the ratio 

Re = ^ (1) 

to quantify this. Thus ii Re ^ 1, the decay is sequential, 
where as if > 1, then the decay is democratic. Large 
values of F/ and thus of Re imply a short-lived interme- 
diate state and thus a difficulty in distinguishing the two 
decay steps. Another way to quantify this is to estimate 
the separation of the fragments from the first step at 
average time of the second decay step, i.e., the decay of 
the the intermediate state. This separation is 

where /i is the reduced mass of products produced in the 
first decay step. If cIe is much larger than the typical 
nuclear diameter (~5 fm), then the two decay steps are 
well separated and cannot influence each other through 
nuclear processes. 

Again we emphasize that the quantities Re and cIe 
do not provide sharp distinctions between prompt and 
sequential processes. For intermediate values of these 
quantities, one may consider a decay process that is ba- 
sically sequential, but where final-state interactions be- 
tween the products from the different decay steps are 
still important. If these final-state interactions do not 
washout sequential signals such as the invariant mass 
of the intermediate state deduced from its decay prod- 
ucts, or the angular correlations between the sequential 
decay steps, then one might still claim that the decay 
has a strong sequential character. However detailed pre- 
dictions for the correlations may require a many-body 
calculation. 

The angular correlations in sequential decay are a con- 
sequence that the system passes through an intermediate 
state of well defined spin. Not all potential sequential de- 
cay scenarios have such correlations, for example if the 
intermediate state is J=0, or the orbital angular momen- 
tum removed in either step is zero. Thus the search for 
these correlations cannot be performed for all cases. 

An example of a truly sequential process is the decay 
of the ground state of to the p+2a exit channel. The 
decay begins by a proton emission to the ground state 
of ®Be which has a very small width (F/=6.8 eV). Here, 
i?g=3.6xl0^^ and ^£=6.1x10^ fm which indicates that 
there are no significant interactions between the proton 
from the first step and the a particles from the second 
step. 

The nucleus ^Bcg s. which decays to the 2p+a channel 
provides an example of democratic decay. The possible 
intermediate state ^i^ig.s. is very wide (Fb^I.23 MeV), 
but most of the strength associated with this state is en- 
ergetically inaccessible in ^Beg.s. decay (see Fig. [3]). Se- 
quential decay is only possible through the low-energy 
tail of this resonance. Thus ^Be^.s. is thus almost a 
Goldansky-type decay. In a sequential scenario, the mean 



energy released in the first step can be determined in an 
R-matrix approximation (Sec. IIIip as i?^=0.64 MeV and 
this gives i?£;=1.92 and dE=G-5 fm, which is not in the 
realm of a truly sequential decay. Experimental studies 
show no evidence of the angular correlations expected 
in a sequential scenario Q (see also Sec. IIV Ap and a 
good description of the decay correlations, the decay en- 
ergy and its width can be obtained with a 3-body cluster 
model 0. 

In this work we report on an experimental investigation 
of particle-unstable states in Be, B, and C isotopes all of 
which undergo disintegration into 3 or more final frag- 
ments. In each case, we have measured correlations be- 
tween these fragments to deduce the nature of the decay 
and to determine to what extent sequential and prompt 
many-body processes are involved. The states reported 
on were all created as projectile-like fragments after the 
interaction of a secondary E/A=70 MeV beam with 
a Be target. One case of particular interest is the Ap+a 
decay of ®Cg.s. As ^Be^.g. is a possible intermediate state 
in the decay of this level, we have obtained improved data 
for the ^Be system using a E/A=70 MeV ''Be secondary 
beam to help in the interpretation of the *C data. Some 
of the ®C results have already been published in Ref. 
The details of the experiment are discussed in Sec. [Hi Re- 
sults for each of the examined states are given in Sec. IIVI 
and the conclusions of this work are presented in Sec. IVl 



II. EXPERIMENTAL METHOD 

A primary beam of E/A—150-MeV ^^O was extracted 
from the Coupled Cyclotron Facility at the National Su- 
perconducting Cyclotron Laboratory at Michigan State 
University with an intensity of 125 pnA. This beam 
bombarded a ^Be target producing i?/^=70.0-MeV ^C 
and ^Be projectile-fragmentation products which were se- 
lected by the A1900 separator with a momentum accep- 
tance of ±0.5%. The ^C secondary beam had an intensity 
of 1.6x 10^ s~^ with a purity of ~65% with the main con- 
taminant being ^Li. The ^Be beam had an intensity of 
4x10^ s^i with a purity of ~90%. 

The two secondary beams impinged on a 1-mm-thick 
target of ^Be. Charged particles produced in the reac- 
tions with this target were detected in the HiRA array 
0. For this experiment, the array consisted of 14 E- 
AE [Si-CsI(Tl)] telescopes located at a distance 90 cm 
downstream from the target. The angular coverage of 
the array is illustrated in Fig. [T] and it subtended po- 
lar angles from 1.4° to 13°. Each telescope consisted of a 
1.5-mm thick, double-sided Si strip AE detector followed 
by a 4-cm thick, CsI(Tl) E detector. The AE detectors 
are 6.4 cmx6.4 cm in area with each of the faces divided 
into 32 strips. Each E detector consisted of four sepa- 
rate CsI(Tl) elements each spanning a quadrant of the 
preceding Si detector. Signals produced in the 896 Si 
strips were processed with the HINP16C chip electron- 
ics [q. For the ^C beam, the time of fiight measured 
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FIG. 1. The angular coverage of the HiRA array used in this 
experiment. 



between a thin scintillator foil in the A1900 extended fo- 
cal plane and the HiRA trigger was used in conjunction 
with the energy loss in this foil to reject most of the beam 
contaminants. 

The energy calibration of the Si detectors was obtained 
with a ^^*Th a-particle source. The particle-dependent 
energy calibrations of the Csl detectors were achieved 
with cocktail beams selected with the A1900 separa- 
tor. These include protons {E/A=60 and 80 MeV), a- 
particles {E/A^60 and 80 MeV), deuterons {E/A=21 
and 59 MeV), tritons (£;/yl=27 and 36 MeV), ^He 
{E/A=37 and 103 MeV) and ^Li {E/A^60 and 80 MeV) 
beams. 

The HiRA telescopes have excellent isotope separation 
for the all of the light fragments of interest in this work 
and multiple fragments within a single telescope can be 
identified The gains for the Si amplifiers were setup 
for the detection of hydrogen and helium isotopes. Some 
information on lithium isotopes is used in this work, but 
lower-energy Li fragments saturated the Si shaping am- 
plifiers and could not be identified (see later). 



III. SIMULATIONS 

The effect of the detector acceptance and resolution is 
quite important when extracting decay widths, branch- 
ing ratios and correlations between the fragments. These 
effects have been extensively studied using Monte Carlo 
simulations that include the angular acceptance and the 
angular and energy resolutions of the detectors, the en- 
ergy loss [l3 and small- angle scattering [TJ] of the frag- 
ments as they leave the target, and the beam-spot size. 
These simulations have proven quite reliable in past ex- 



periments with HiRA [9',(T1. 

A number of different types of simulations are de- 
scribed in this work including both sequential and 
prompt decay processes. However, the predicted reso- 
lution of the excitation energy of a level was found to be 
mostly insensitive to the details of the decay process, but 
is largely determined from the detector properties. 

In simulating sequential 3-body decays, where the 
width of the intermediate state is significant, we have 
used the R-matrix formalism [ij, to predict the dis- 
tribution of Et, the total kinetic energy released, and 
Ex=Et-EI, the kinetic energy released in the second 
step; 



NiET,E,) 



TiIEt^E,) 



[Et -Qi- AUEtW + (l/4)rL(i?T) 



where 



TtotiEr) = J Ti{ET,E^)dE,, 



oo 

TtotiEr) = j Ai{ET,E,)dE^, 



Ti{Et,E,) = 2-ilPu{ET - E^)p{E,), 



(3) 



(4) 



(5) 



(6) 



Ai{Et,E,) = --f^,[Su{ET ~ E,) - Su{Qi - E^ME,), 

(7) 

(8) 



p(E^) = c 



r2(s.) 



[E^ - Q2 - A^iE^W + {l/4)TliE,)' 

OO 

J p{E^)dE^ = 1, 



V2{E^) = 2jiP2i{E,), 
ME.) - -"fi[S2e{E,) - S2i{Q2)l 



ka 



Si{E) 



_ Fl,{ka)Fi{ka) + G'f,{ka)Ge{ka) 
" Fi{kay + Ge{kay '' 



lA5im{AY^ + A. 



1/3n 



(9) 

(10) 
(11) 
(12) 

(13) 
(14) 



F and G are regular and irregulat Coulomb wavefunc- 
tions, k is the wave number, c is a normalization con- 
stant determined from Eq. -ff and 7I are the reduced 
widths associated with the first and second decay steps, 
and Qi and Q2 are the centroids associated with Et and 
Ex, respectively. Angular correlations in sequential de- 
cay are calculated from Refs. [H, [l^. 
The reduced width can be expressed as 



fi2 



(15) 



where S is the spectroscopic factor and 9^^, the single- 
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particle dimensionless reduced width, is 
o2 a '^^(a) 



2 J^u^{r)dr' 



(16) 



Here u{r) jr is the single-particle radial wave function cal- 
culated with a Coulomb plus Wood-Saxon potential with 
standard parameters for radii and diffuseness (ro=1.25 
fm, rc=l-3 fm, and a=0.65 fm) and the depth adjusted 
to fit the resonance energy. 



IV. RESULTS 

The deduced properties of the levels investigated in 
this study are summarized in Table H] including their cen- 
troids, widths, and decay modes. More detailed discus- 
sion for each case are contained in the rest of this section 
including comparison with the evaluated quantities from 
the ENSDF database Excitation energies are de- 

termined from the invariant mass method, i.e., the total 
kinetic energy of the fragments in their center-of-mass 
reference frame minus the decay Q value. 



C ground state 



i2 200 



(a) 4p+a 




The ground state of ®C is unstable to disintegration 
into four protons and an a particle. The distribution of 
®C excitation energy reconstructed from the Ap+a chan- 
nel following neutron knockout of beam particles is 
displayed in Fig. [D^a). The excitation energy was de- 
termined assuming the mass excess of 35.094 MeV 
obtained from the 2003 mass evaluation of Ref. [l^ ■ The 
mass excess is needed to calculate the decay Q value. 
The excitation-energy spectrum shows a strong peak lo- 
cated near zero excitation energy indicating that we have 
populated the ground state of ®C. The long-dashed curve 
shows the simulated shape, including detector response, 
based on the evaluated mass excess and the listed *C 
decay width of 230 keV. This curve also contains a back- 
ground contribution which is indicated by short-dashed 
curve. The simulated shape is both wider and shifted up 
in energy compared to the experimental peak suggesting 
that the evaluated mass excess and width used in the 
simulations are incorrect. 

In order gain a better understanding of the magnitude 
of the experimental uncertainties associated with this 
measurement, we have investigated the line shape of the 
^Be ground state formed in neutron knockout reactions 
from the ^Be beam. In the and ''Be neutron knock- 
out reactions, the beam velocities, target thickness, and 
detector apparatus were identical and thus systematic 
errors are expected to be the same. The reconstructed 
^Be excitation-energy spectrum from the detected 2p+a 
events is displayed in Fig.[2lb). This spectrum contains a 
peak at zero energy associated with the ground state and 
a wider peak at ~1.7 MeV associated with the first ex- 
cited state. For the ground-state peak, the solid curve in 



FIG. 2. (Color online) Excitation-energy spectra for (a) *C 
fragments produce in neutron knockout from the beam and 
(b) ®Be fragments produced in neutron knockout from the ^Be 
beam. In both cases, the curves show simulated ground-state 
distributions (see text). The solid curve in (a) is a fit to the 
data with a background indicated by the dotted curve. The 
dashed curve is the simulated response based on the tabulated 
peak location and width. In dotted and dashed curves in (b) 
are the two possible background distribution considered in the 
text. 

this figure shows the simulated spectrum, again including 
a background contribution (including a contribution from 
the first excited state) which is indicated by the short- 
dashed curve. The tabulated ®Be mass excess and decay 
width were used in this simulation which reproduces the 
experimental result very well, strongly suggesting that 
the disagreement of the simulation with the data for *C is 
not an experimental artifact. The experimental FWHM 
of the ^Be ground-state peak is 214±10 keV, significantly 
greater than its intrinsic width of r=92±6 keV. Thus the 
agreement between the ^Beg.s. data and simulation indi- 
cates that we can correctly account for the experimental 
resolution. 

The solid curve in Fig. [2ja) shows the results of a simu- 
lation where both the *C mass excess and its decay width 
were adjusted to best reproduce the experimental peak. 
The fitted mass excess is 35.030±0.030 MeV and the fit- 
ted width is 130±50 keV. 

Previous measurements of the mass excess and decay 
width were made in Refs. [T9l - [2l| using transfer reactions. 
Of these, only the 1976 work of Tribble et al. |21| has 
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TABLE I. The excitation energy, widtii, spin and isospin of states for which new information is determined. 



nucleus 


decay 


branching 
ratio 


exit 
channel 


E* 
[MeV] 


r 

[keV] 




T 


TT^ 


H-r^ 

P + ^BSg.^. 


81±10% 


3p+Q 


0.0 


801±20 




3/2 


^B 


P+ Be4.57 




p+'^He+a 


5.93±0.02 


850±260 


3, 4 


1 


«B 


P+''Be6.73 




p+^He+a 


8.15±0.20 


950±320 




1 


®B 


2p+^Li2.18 

2p+d+a 


~97.5% 
~1.0% 
~1.5% 


2p+^Li 
2p+rf+a 


10.619±0.009^ 


<60^ 


0+ 


2 


R-i-i 

'^Be 


p+ Ll4.63 




p+t+a 


22.96±.02 


680±146 


3, 4 

0+ 


1 


«C 


2p+''Beg... 


100% 


4p+a 


0.0 


130±50 


2 


^B 


a+'^Li 
p+ Bea.os 


97.2±0.5% 
2.8±0.5% 


p+2a 
p+2a 


11.70±0.02 


880±80 


3/2-,5/2+,7/2- 


1/2 


^B 


p+^BeT=o+i 




p+2a 


16.99±0.03 


22±5^ 


1/2- 


3/2 


''B 


p+*Bei9.069 




2p+'^U 


20.64±0.10^ 


450±25'o 






10c 


P+^B2.345 




2p+2a 


8.54±0.02~ 


<200 






10(j 


2p+2a 
a+'^Be 

P+^B2.34 


48% 
35% 
17% 


2p+2a 
2p+2a 
2p+2a 


9.69 


490 




1 








2p+2a 


10.48±0.2 


<200 






10^ 






2p+2a 


11.44±0.2 


<200 







^ Assumed decay is to the ground state of ^Li. 
value from tabulations 



more than a handful of ®C events. The reaction stud- 
ied was ^^C('*He,^He)^C and the *He fragments were de- 
tected. The extracted mass excess was 35.10±0.03 and 
the decay width was 230±50 keV and 183±56 keV as- 
suming Breit-Wigner and Gaussian intrinsic lines shapes, 
respectively. As in the present work, the experimental 
resolution was significant. 

The *C ground state is quite unusual in having a 
five-particle exit channel. It is of interest to determine 
whether this state disintegrates directly into five pieces 
or proceeds in a sequential process of two or more steps 
through intermediate states. Figure [3] shows the levels of 
possible interest in the decay. Of the possible intermedi- 
ate states, the ^Be ground state is the narrowest (longest 
lived) and easiest to detect experimentally. In each de- 
tected Ap+a event, one can find 6 ways to reconstruct a 
^Be decay as a 2p+a subset. For each of these ways, a 
^Be excitation energy is determined. We use gate Gl in 
Fig. [2][a) to select out the ground-state events. The re- 
sulting excitation-energy spectrum of all these potential 
^Be candidates is displayed in Fig.l^a) as data points. As 
there can be at most only one ^Be ground-state fragment 
per event, at least 5/6 th of the spectrum corresponds to 
incorrectly reconstructed ^Be fragments and these form 
a background. The experimental spectrum clearly dis- 
plays a peak centered near zero energy associated with 
the ^Be ground state on a background of approximately 
Gaussian shape. 

To determine what fraction of the ^G events decay 
through ^Beg,s., we have fitted the spectrum in Fig. Ilja) 
with two components; a ^Be^.s. component with shape 
taken from the experimental peak in Fig. [2][b) obtained 
with the ^Be beam and a background component taken to 
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FIG. 3. Levels of interest in the decay of the *C and ^B 
ground states. 



be Gaussian in the immediate vicinity of the ^Bg.s. peak. 
In order to use the experimental ^Beg.s. line shape from 
Fig-EIb), one has to first subtract the small background 
under this peak. The short and long-dashed curves in 
this figure show two possible backgrounds which were 
considered. The short-dashed curve is associated with 
the solid curve obtained using a symmetric Breit-Wigner 
line shape. The use of the long-dashed background curve 
would imply a slightly asymmetric ^Be^ line shape 
which has an enhanced higher-energy tail. A small degree 
of asymmetry is not unreasonable. 

The resulting fit (using the original ^Be background) is 
shown by the solid curve in Fig. IH^a) where the Gaussian 
background is indicated by the dashed curve. Taking into 
account both possible ^Be backgrounds, these fits imply 
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FIG. 4. (Color online) The spectrum of all possible recon- 
structed ^Be excitation energies for events in the ground 
states of (a) *C and (b) '^B. The data points are the ex- 
perimental data, while the solid curves show fits to the ^Be 
ground-state region with the dashed curves showing the fitted 
background. 



that, on average, 
ated in each 



g.s. 



1.01±0.05 Beg s fragments are cre- 
decay, i.e. essentially all decays pass 



through ^Bcg.s.-This number was corrected for the small 
background under the ^C^.^. peak itself [dotted curve in 
Fig. Elja)], using the the region Bl to estimate the ^Be 
probability associated with this background. 

In many ways, the decay of ^Cg,s to 2p+^Beg,s. is 
similar to the decay of ^Beg s. to 2p+a. In each case, 
both the initial and final states are and the pos- 

sible intermediate states CBg,s. and ^Lig.s) are J^=Z/2~ . 
These intermediates states are both wide, but only their 
low-energy tails are energetically accessible. Both decays 
are democratic; i?B=1.2 and 1.9 and (i£;=8.8 and 6.5 fm 
for *C and ^Be, respectively. Thus the disintegration of 
^Gg,s. is novel having two sequential steps of democratic 
two-proton decay. To investigate this possibility, it is use- 
ful to separate out experimentally the protons from the 
first and second decay steps. However, it is not possible 
to do this with 100% certainty for all decays. 

For a fraction of the events it is possible to make this 
separation with reasonable precision. To this end, we 
have selected events where one, and only one, of the 6 
possible reconstructed ^Be fragments has an excitation 
energy associated with the ground-state peak indicated 



by the region G2 in Fig. |4l^a). For the events which sur- 
vive this criteria, the identities of the ^Bcg.^. decay prod- 
ucts are obvious. 

Some events from the tails of the ^Bcg.s. peak lie out- 
side the G2 region and will be rejected. Also, and more 
importantly one can have the situation that the correct 
^Be fragment is in the high or low-energy tail outside 
of the gate and an incorrectly identified ^Be fragment 
lies within the G2 gate. Such an event will not be re- 
jected, but the identities of the protons from the first 
and second steps will be incorrect. Monte Carlo simula- 
tions were used to optimize the width of the G2 region so 
as to obtain the best compromise between event-rejection 
and mis-identification. 

With the chosen gate width we select only 30% of the 
^Gg,s. events and, of these, we estimate that a selected 
event is incorrectly identified 30% of the time. While this 
is significant, these misidentifications involve cases were 
the two protons assigned to one of the decay steps are 
actually from the two separate steps. As the two steps 
are separated by a ^Bcg s. fragment with J'^^O^, there is 
no angular correlations between the protons from the two 
steps and these misidentified events will give rise only to 
smooth backgrounds in the correlation plots. 

The energy and angular correlations between the par- 
ticles produced in two-proton decay can be described by 
the hyperspherical Jacobi vectors X and Y and their con- 
jugate momenta k^^ and kj,. There are two independent 
ways of defining the coordinates which are referred to as 
the "T" and "Y" systems. These are illustrated in Fig. [5] 
where the core (fragment 3 in the "T" system or frag- 
ment 2 in the "Y" system) is the ^Be^.^. fragment in the 
first 2p decay and the a particle in the second 2p decay. 
In terms of the position vectors r^, momentum vectors 
ki and masses mi 2, and 3) the Jacobi coordinates 

are 



X = ri - r2 , 

miVi + m2Y2 



mi - 
m2ki 



m2 

- mik2 



mi -I- 7712 

m3(ki -I- k2) - (mi + m2)k3 

mi + m2 + m3 



(17a) 
(17b) 

(17c) 

(17d) 



Of the six degrees of freedom required to define the k^; 
and ky distributions, three describe the Euler rotation of 
the decay plane and one is constrained from energy con- 
servation. Thus the complete correlation information can 
be described by two variables which we take as Ex/Et 
and Ok, where E^ is the energy associated with the X 
coordinate. 



Ex — 



(mi -I- m2)kl 



(18) 



2mim2 

Et is the total three-body energy and 6k is the angle 
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FIG. 5. Independent "T" and "Y" Jacobi systems for the 
core+N+N three-body system in coordinate and momentum 
spaces. 



between the Jacobi momenta, 



For each event, there are two ways of labeling the two 
protons and thus there are two possible values of the 
[Ex/ ET,cos{9k)] coordinate, both of which are used to 
increment the correlation histograms. For the "T" sys- 
tem, this produces a symmetrization of the angular dis- 
tributions about cos(0fc)=O. 

Distributions constructed in the two Jacobi systems 
are just different representations of the same physical pic- 
ture. Each Jacobi system can reveal different aspects of 
the correlations. The Jacobi "Y" system is particularly 
useful if the two protons are emitted sequentially through 
an intermediate state. In this case, is the kinetic en- 
ergy released in the second step and 6k is the angle be- 
tween the two decay axes. The Jacobi "T" system can 
have a similar interpretation for a diproton decay with 
an intermediate ^He state. 

The projected correlations in both the "T" and "Y" 
systems are plotted in Fig. [5] for the second decay step of 
^Cg,s., i.e., the 2p decay of the ^Be^.s. intermediate state. 
These can be compared to results from the ''Be beam in 
Fig. [3 where ^Be^.s. is formed more directly following 
a neutron-knockout reaction. In this case, the statisti- 
cal uncertainties are significantly smaller and there are 
no misidentified events. A consistency between the pro- 
jected correlations in Fig. [HI and [7] is necessary if the ^Cg.s. 
level disintegrates through ^Be^.s. and indeed the corre- 
sponding projected correlations are quite similar in the 
two figures. 

In order to make more detailed comparisons of the ex- 
perimental ^Be correlations in Figs. [H] and [71 we have 
made extensive Monte-Carlo simulations including the 
effect of the detector response and the event selection re- 
quirements. The results from the ^Be beam in Fig. [7] are 
compared to simulations where the effect of the detector 




FIG. 6. (Color online) Projected correlations in (a,b) "T" 
and (c,d) "Y" Jacobi systems for the three-body decay of ®Be 
formed in the second step of *Cg.s. decay. Angular correla- 
tions are shown in (b) and (d) and energy correlations are 
shown in (a) and (c). The data points are the experimen- 
tal data, while simulated results are indicated by the curves 
which pass through the data points. The lower curves are es- 
timates of the "backgrounds" of incorrectly identified events. 
The short and long-dashed and solid curves are associated 
with simulations (A), (B) , and (C), respectively (see text). 



resolution and acceptance is considered with no misiden- 
tification or selection bias. The detector effects induce 
only small modifications to the projected correlations in 
this case. The solid curves are predicted distributions 
from the Quantum-Mechanical 3-body cluster model of 
Ref. [H (Calculation P2 in that reference). This model 
reproduces the experimental total energy Et, F, and the 
experimental correlations in Fig. [71 These predicted cor- 
relations will therefore be used in the simulations in the 
disintegration of *C. 

The dashed curves in Fig. [71 are from an R-matrix ap- 
proximation of sequential 2-proton emission from ^Bcg.^. 
through the ^Lig.g. intermediate state as described in 
Sec. mil The angular correlations in this simulation are 
particularly strong and favor decays where the two se- 
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system is taken to be the same as for ^Be^ s. decay. 
Note, as the core mass is different in the ^C^.s. 2- 
proton decay, the Jacobi-"Y" distribution wih not 
be exactly the same as for ^Be^.s. decay. 

(C) The decay is treated in the R-matrix approxima- 
tion as a sequential 2-proton decay through the ''B 
intermediate state in a similar manner as just de- 
scribed for the sequential-2 proton decay of ^Bcg 
The R-matrix parameters for ^Bg.s. are taken from 
Sec. UVBl 

The short-dashed, long-dashed, and the solid curves 
in Fig. |6] are the results for ^Bcg.s. decay in the second 
step oi^Cg,s. disintegration for simulations (A), (B), and 
(C), respectively. The curves, which pass through the 
data points in each panel, include the detector response, 
the selection bias and the contribution from the incor- 
rectly identified events. The latter contributions are in- 
dicated by the lower curves in each panel. In all cases, 
the three curves from the different simulations almost 
overlap. Thus the simulated results for ^Be^.^. decay 
show very little sensitively to the decay correlations in 
the first ^C^2p+^Beg,s. step. Our assertion that the 
distributions for the misidentified events show no strong 
structure is borne out and the simulations for all se- 
lected events reproduce the experimental data. There- 
fore within the statistical errors on the experimental data 
and uncertainty associated with extracting the data, that 
correlations in the second step of disintegration is con- 



sistent with "Be^.s. decay. 

For the first step {2p+^Beg,s.) of ^C^.^. decay, the pro- 
FIG. 7. As for Fig.H but now for ^Be, , events produced via Jected E^/Et and cos(6'fc) distributions for both the Ja- 

cobi "T" and "Y" systems are plotted in Fig. E The 
curves in each panel show the simulated distributions and 
backgrounds from the three simulations discussed above. 
Again these backgrounds are similar in all the simula- 
tions and show relatively smooth behaviors. Clearly, 
the phase-space simulation (short-dashed curve) is in- 
consistent with all the experimental projected correla- 
tions. Also the simulations using the theoretical ^Be cor- 
relations (long-dashed curves) does not fit the data in 
Figs. [5]^a) and Hljd) indicating that two-proton correla- 
tions are different in the first and second steps. 

The sequential decay simulations gives the best de- 
scription to all the experimental distributions, however, 
it is clearly not perfect. The predicted angular correla- 
tions indicated by the solid curve in Fig.[8l^d) show strong 
enhancements near cos(0/c)=±l. These angular corre- 
lations also cause the large structure predicted for the 
Jacobi "T" energy distribution of Fig. |5{a). The experi- 
mental data also show an enhancement near cos(0fe)=-l, 
sometimes called the diproton region, however, the cor- 
responding enhancement near cos(0)=l is not evident. 
In the sequential scenario, cos(0)=l corresponds to the 
®Be core from the second step being directed towards 
the proton from the first step. Coulomb final-state in- 
teractions may deflect the trajectories and thus suppress 
the yield near cos(0)=l. Alternatively the enhancement 



neutron knockout from the ^Be beam. The solid and dashed 
curves show predictions of the 3-body cluster model and a 
sequential-decay calculation, respectively. 



quential decay axes are collinear, i.e. cos{9k) = ±1 in 
the Jacobi "Y" system [see Fig. (TJd)]. This angular cor- 
relation is also responsible for the predicted broad mini- 
mum in the energy distribution in the Jacobi "T" system 
[Fig. [Tl^a)]. However, these features associated with the 
predicted angular correlation are lacking in the experi- 
mental data. This indicates that the decay of ^Be^.s. does 
not pass through an intermediate state of well defined an- 
gular momentum (J'^=3/2~). Similar conclusions have 



been made in other studies of ^Be 



g.s. taiL£iaj- 

Now let us return to "^Be^.s. decay in the process 
of ^Cg,s. disintegration. In order to include the bias 
from the event selection and the contribution from 
misidentified events, it is necessary to simulate the first 
step of ^Cg.s. disintegration. For this first decay step, 
^C— !'2p+^Bcg.s., we have considered three quite differ- 
ent sets of simulated [Ex/Et, cos(6'fe)] correlations: 

(A) The momenta of the three fragments are chosen ac- 
cording to available phase-space volume. 

(B) The [E.x/Et, cos{0k)] correlations for the Jacobi "T" 
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near cos(6'fc)=-l may not be related to a sequential angu- 
lar correlation but rather might be a signature of some 
"diproton" character of the decay similar to that ob- 
served for the 6.57-MeV state in 

The total decay width in the sequential scenario can 
be estimated as [3, [T^ 



1 + 7?/; 



oo dSit(ET- 
dEx 



' Ex=Q 



p{Ex)dE, 



(20) 



The p+^Bg.s. spectroscopic factor was calculated with 



CKI Hamiltonian 
the center of ma 



124] as S'=3.20 (8/7) 
jS correction. With 



where (8/7) is 
6)2p=0.44 from 
Eq. ([H]) we estimate 7^=4.28 MeV and thus r°=12 keV 
which is significantly less than the experimental value of 
130±50 keV and thus is clear that this sequential de- 
cay scenario cannot provide a total description of the 
2-proton decay of ^Cg.s.. 

Given the enhanced "diproton" character of the decay 
it is useful to consider the diproton decay model which 
can be treated in the R-matrix formulation of Barker (25j 
together with the diproton cluster decay model developed 
in Ref. [1^. For the wavefunctions we use the p-shell 
basis with the CKI Hamiltonian '24'] . 

The diproton decay spectroscopic factor is given as a 
product of the three terms as in Eq. (7) of [1^. For 
the ®C decay these are (8/6)^ = 1.78 (the center of mass 
correction), G^{p)=l/2 (the cluster overlap factor) and 
C=1.21 (the p-shell spectroscopic factor for L—S=Q), to 
give 562=1-08. The related single-particle dimensionless 
reduced width of 6*^^=1.00 is also taken from [2^. The 
diproton decay width obtained for the experimental de- 
cay energy of Qi=2.077 MeV is r=88 keV. If we use 
the older Qi value of 2.147 MeV based on the tabulated 
masses, the width would be 100 keV. We have assumed a 
narrow ^Be final state. The change coming from a folding 
with its actual width is small. Our spectroscopic factor 
Se2 is a factor of five larger than given by Barker fol- 
lowing Eq. (13) of Ref. [2g]. We do not know why, but 
assume that there was an error in Barker's calculation. 

The calculated diproton decay width of r=88 keV is 
consistent with the experimental value of 130±50 keV, 
supporting the interpretation of the enhancement at low 
Ex/Et in the Jacobi "T" system [Fig.|S{a)] as a diproton 
feature. However, the diproton model does not provide 
a complete description of the experimental correlations. 
For example in the diproton model, the Jacobi "T" an- 
gular distribution is just the distribution of the diproton 
decay angle, which should be isotropic as the diproton 
has J^=0'^. However the experimental distribution in 
Fig. ^ih) is not uniform and thus inconsistent with this 
notion. 

We have also considered the diproton model for ^Bcg s. 
decay. For the spectroscopic factor we obtain (6/4) ^ = 
2.25 (the center of mass correction), G'^{p)=\/2 (the clus- 
ter overlap factor) and C=1.00 (the p-shell spectroscopic 
factor for L=S'=0), to give 542=1.12. With the single- 
particle dimensionless reduced width of 9^^^=1.121 [26| . 




FIG. 8. As for Fig. El but now the correlations are associated 
with the first step of *Cg.s. decay forming 2p+®Beg.s.. 



the diproton decay width obtained for (31=1.371 MeV 
is F=98 keV again consistent with the tabulated value 
of 92±6 keV. However in this case the diproton nature 
of the experimental correlations (Fig. [7]) is even less ap- 
parent. Clearly a complete understanding of ^Be^.s. and 
^Gg.s. requires a full three-body decay calculation. 



B. 



levels 



The ^B excitation-energy spectrum determined from 
all detected ip+a events is displayed in Fig. [SKa). The 
peak near zero excitation energy indicates that the 
ground state of ''B was produced in the reaction. There is 
a significant contamination in this spectrum from ^Cg.s. 
disintegrations where only 3 of the 4 protons were de- 
tected. To determine the contribution from this pro- 
cess, we have taken detected ^p+a events associated 
with ®Cg.s. and generated pseudo 'ip-\-a events by alter- 
natively removing one of the protons and determined the 
^B excitation energy. The distribution of these excita- 
tion energies is labeled as ^Cg.s. in Fig.lHl^a) and indicates 
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FIG. 9. (Color online) (a) The excitation-energy spectra for 
events determined from all detected 3p+Qf events and from 
^Cg.s. 4p+a events where one of the protons was ignored, (b) 
The data points are the excitation-energy spectrum after 
the *C contamination was subtracted. The solid curve shows a 
fit to this data assuming an R-matrix line shape and including 
the detector response. The long-dashed curves indicate the 
background obtained from the fit and the raw R-matrix line 
shape is indicated by the dotted curve. 



that the incompletely detected *Cg.s. events lead to an 
enhancement of the low-energy side of the '^Bg.s. peak. 
The Monte Carlo simulations were used to normalize the 
^Cg,s. distribution to the expected number of incomplete 
detected events. The data points in Fig. Mj^) show the 
results for the excitation-energy spectrum after this 
^Cg,s. contamination was subtracted. 

The subtracted distribution peaks at E* ^-100 keV 
rather than zero, suggesting that the evaluated mass ex- 
cess from Ref. [l^ used to calculate the Q value for 
this figure is incorrect. This evaluated mass excess of 
27.94±0.10 MeV was derived from an average of two 
measurements; 27.80±0.10 MeV from the ^Li(7r+,7r")^B 
reaction and 27.94±0.10 MeV from the i°B(3He,6He)'^B 
reaction p7l |. 

As this is a rather wide level near threshold, the use 
of a Breit-Wigner line shape r nay not be appropriate. 
We have used R-matrix theory [iJl assuming p-t-^Be^.^. 
as the only open decay channel. The line shape is 
given by Eq. ([5]) where now is the energy above the 
p+^Beg,s. threshold, and Q2 is the resonance energy. Us- 



ing J''=3/2- listed in the ENSDF database and £=1, the 
experimental spectrum was fit by adjusting E2 and 7I 
and other parameters defining a smooth background un- 
der the peak. The effects of the detector acceptance and 
resolution was taken into account via the Monte Carlo 
simulations. The resulting fit is displayed as the solid 
curve in Fig. and the fitted background is shown as 
the dashed curve. The dotted curve shows the fitted peak 
shape without the effects of the detector acceptance and 
resolution. 

The fitted resonance energy is Q2=2.013±0.025 McV 
which implies a mass excess of 27.677±0.025 MeV. The 
error includes the effect in a 20% uncertainty in the mag- 
nitude of the subtracted ^C contamination. This mass 
excess is consistent with the previous measurement from 
the ^Li(7r+,7r^)^B reaction, but is inconsistent with that 
from the ^°B(^He,^He)^B reaction. The fitted reduced 
width is 7|=1.32±0.02 MeV and thus the level widthjH], 



nQ2 



1 



^2 dSe. 
'2 dE^ 



(21) 



is 0.80±0.02 MeV. The evaluated width in the 
ENDSF database is 1.4±0.2 MeV derived from the 
^°B('^He,^He)''B reaction 27[. Again our results are in- 



consistent with this measurement. 

In the above analysis of this level, we assumed ^Bg ^ — > 
p-t-^Bcg.s. sequential decay. We can check this by looking 
for a ^Beg fragment as described in Sec. IIV Al Here, 
there are only three ways to reconstruct a ^Be fragment 
from 3 protons and an a particle. The distribution of 
all three possible ^Be excitation energies is shown in 
Fig. Ub) for the gate G3 in Fig. M^} on the ^B exci- 
tation energy. This gate does not cover the whole of the 
^Bg s peak width, but avoids the lower energies where the 
*C contamination is present. The ^Be excitation-energy 
spectrum displays a prominent peak at zero excitation 
energy indicating that ^Be^.s. fragments were produced 
in the decay. The solid curve shows a fit to the data us- 
ing the experimental ^Be^.s. line shape [Fig. [5l^b)] and a 
smooth Gaussian background (dashed curve) under this 
peak. From this fit we determine that, on average, there 
is a 54±6% probability of finding a ^Be^.^. fragment in 
the G3 gate. However it is clear from Fig. [IJb) that 
there is a very significant background under this peak. 
Using gate B2 in Fig. EJb) we estimate that the back- 



ground has a 19±2% probability of containing a ^Be 
fragment. Using this, we find that 81±10% of the ''B 



g.s. 
g.s. 

events decay via p-t-^'Be^.g.. 

Shell-Model calculations for ''Bg.s. with the CKI 
Hamiltonian [IJl give the spectroscopic factor for the 
p+^Be configuration as S'=0.59(7/6) = 0.688. The 
p-|-^Bei.67(2''") configuration is predicted to be '^3 times 
stronger but decay to this channel is suppressed due to 
its significantly smaller barrier penetration factor. With 
6*3^=0.688 from Eq. (fTS]) . this gives us a predicted re- 
duced width of 1.42 MeV [Eq. ^^]. Taking into account 
the 81% branching ratio, we estimate the experimental 
value to be 7^=1.07±0.15 MeV which is quite similar to 
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FIG. 10. (Color online) Level scheme showing the levels of *B 
discussed in this work and their decay paths. 



the predicted value confirming that the p+^Tig,s. is not 
the strongest configuration in ^Bg . 

Although it is possible to consider a proton decay 
through the 1.670-MeV (J''=2+) first-excited state of 
^Be (see Fig. [3]), this ^Be intermediate state is sufficiently 
wide that the disintegration is democratic. The average 
energy released in the first step is £'^=343 keV (deter- 
mined from the differences in centroids of the two levels) 
which is very small compared to the width of the inter- 
mediate state, r=1.16 MeV. Alternatively, the distance 
traveled by the proton on average before the ®Be frag- 
ment decays is 2.0 fm which is small compared to the 
nuclear diameter. Thus this decay strength is probably 
best described as a 4-body decay. 



C. levels 

The ^B excited states studied in this work and their 
observed decay paths are illustrated in the level diagram 
of Fig. [TOl We observed excited states in three detected 
exit channels; p+^He+a, 2p+^Li, and 2p+d+a. The 
p+^Be channel, which is probably the most important 
exit channel for many states, was not accessible in this 
work as Be fragments saturate the Si shaping amplifiers. 



1. p+^He+a exit channel 

The ^B excitation-energy spectrum derived from all 
detected p+^Yie+a events, shown in Fig. [TlT a). displays 
wide peaks at 5.93 and 8.15 MeV. The ^Be excitation- 
energy distribution obtained from the ■^He-l-a pairs asso- 
ciated with each of these events is shown in Fig. lTlT b). A 
peak associated with the 4.57-MeV 7/2~ level is promi- 
nent. The total ®B excitation spectra is subdivided in 
Fig. ITlT a) into those in coincidence with the ^Be peak 
[gate GA in Fig. ITlT b)] and those that are not. This sub- 
division separates the two ^B peaks. The 5.93-MeV peak 
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FIG. 11. (Color online) Excitation-energy spectra determined 
from p-|-^He-|-a events, (a) The *B excitation spectrum from 
all p+"^IIe-|-Qf events and its decomposition into those asso- 
ciated with the 4.57-MeV ^Be intermediate state and those 
not. (b) The ^Be excitation spectrum obtained from the '^He- 
a pairs of all p-|-'^He-|-a events, (c) Background-subtracted 
^Li excitation-energy spectrum associated with the 8.15-MeV 
state in *B. Vertical arrows locate known levels with their 
widths indicated by the error bars. The gates used in this 
work are also indicated by the dashed lines. 



is associated with proton decay to the 4.57-MeV state of 
''Be and 8.15-MeV state does not decay by this path. Let 
us first concentrate on the 5.93-MeV peak. 

As there are only three fragments in the exit channel, 
the velocity vectors of the decay products are located 
in a plane in the *B center-of-mass frame. To display 
the correlations between the decay fragments, we have 
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FIG. 12. (Color online) Distributions of velocities of the three 
fragments for the p+'^He+a exit channel associated with the 
5.93-MeV *B state. Events were selected using the G4 and 
G6 gates in Fig. 1111 The distributions are projected onto the 
plane of the decay, with all "^He and all a-particles locations 
forced to be approximately coincident. 
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projected these vectors onto that plane, and then within 
this plane, rotated all the vectors such that the locations 
of all the Qf-particle velocities approximately coincide and 
at the same time the locations of all the ^He fragments 
velocities approximately coincide. This was achieved by 
requiring that the relative ^He-a velocity was parallel to 
the Vz axis. A contour plot of the resulting distribution of 
p, ^He, and a velocities is displayed in Fig. [T^l The ■^He- 
a separation is constant as these fragments come from 
the decay of the 4.57-MeV ^Be state. The protons lie 
approximately on an arc centered at the origin (^B center 
of mass) as indicated by the dashed curve. Clearly the 
protons were emitted in the first step and the magnitude 
of their velocity is independent of their emission direction 
and the subsequent decay axis of the ^Be intermediate 
state. 

The decay of this state is expected to be sequential 
with i?s=0.11 and dE=(yO fm, and thus the two decay 
steps are largely independent, apart from consideration 
of conservation laws. Angular momentum conservation 
gives rise to the angular correlations. Let us consider the 
relative angle 9i2 between the two decay steps which is 
related to the angle 6k in the Jacobi "Y" system. A O12—O 
corresponds to an event where the proton in the first step 
and the a in the second step are emitted in the same di- 
rection. The distribution of 612 is shown as the data 
points in both panels of Fig. 1131 The experimental dis- 
tribution is symmetric about cos(0i2)=O consistent with 
the expectation for a sequential decay scenario. The ex- 
perimental distribution is clearly not isotropic (flat). The 
predicted correlations in a sequential decay depend of the 



FIG. 13. (Color online) Angular correlations for the 5.93- 
MeV state of *B. Events were selected using both the G4 
and G6 gates in Fig. 1111 The angle 612 in the relative angle 
between the initial proton decay axis and the subsequent ^He- 
a decay axis in the center-of-mass frame. The proton 
and Q particle are emitted in the same direction when 612— 0. 
The solid curves are predictions for pure values of jp, while 
the dashed curves are for mixed values where the degree of 
mixing was chosen to best fit the data. 

initial spin J and the total angular momentum removed 
by the proton jp in the first decay step. The predicted 
correlations do not depend on the parity of the initial 
state for these proton decays, but mixing between jp val- 
ues should also considered, i.e., the angular correlation 
becomes 

w{9i2) = |ai|^ wi(6'i2) -I- I (22 1^^2(6*12) 

+ {aia2+ala2)winter{0i2), (22) 

where wi and W2 are the correlations associated with two 
pure values of jp, ai and 02 are their complex amplitudes 
and Winter is the interference term. 

With this mixing, only two possible J values were con- 
sistent with the experimental distributions. Figure [T3l a) 
shows the results for J=3. The two solid curves show the 
predictions for two pure values of jp; 1/2 and 3/2. These 
predictions have been normalized to the same number 
of events as in the experimental distribution and Monte 
Carlo simulations indicate that the distortions due to the 
angular acceptance and energy resolution of the detectors 
are minimal. Neither of these predictions fit the data at 
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all. However, a mixed solution with relative amplitudes 
of ai=l and a2=0.49±0.14 for the two jp values, respec- 
tively, fits the data very well [dashed curve in Fig. fTSl a)]. 

The other solution was obtained for J=4 and is shown 
in Fig. [TSl' bV Again, pure jp=l/2 and 3/2 predictions 
are shown by the solid curves. In this case the jp—3/2 
predictions fits the data reasonably well. If we allow for 
mixing, the best fit is obtained with relative amplitudes 
of q;2=1 and ai~-0.07±0.10, however, the level of mixing 
is minimal. All told, the angular correlations indicate 
that the 5.91-MeV state is either J=3 or 4. 

For the 8.15 MeV state in ®B, which does not proton 
decay to the 4.57-MeV ^Be state [Fig. E^a)], it is dif- 
ficult to clearly define its decay path due to the larger 
background under this peak and also because there are 
no possible narrow intermediate states through which 
it can decay. Sequential decay through a higher-lying 
^Be state and through the ^Li ground state are possi- 
ble. Apart from the narrow 4.57-MeV state, the ''Be 
excitation spectra in Fig. [TlTb) shows a broad peak at 
-6.5 MeV, close to the location of the £:*=6.73-MeV, 
J-"=5/2-, r=1.2-MeV state in '^Be, indicated by the ar- 
row in the Fig. [TlT b). For a wide intermediate state, 
the full width of the state is not always populated by 
sequential decay as the Coulomb barrier of the first step 
can suppress the higher excitation region. The peak at 
6.5 MeV may be associated with the 6.73-MeV state and 
gating of this peak enhances the 8.15-MeV peak. How- 
ever for all the p+^lle+a events, only 30±6% of the ki- 
netic energy, in the reconstructed *B frame, is associated 
with the proton. Thus the ^He-I-Q; pair accounts for the 
remaining 70±6% and thus the ■^He-a reconstructed ex- 
citation energy is strongly correlated to p+^He+a exci- 
tation energy suggesting that the 6.5-MeV peak in ''Be 
spectra in Fig. [TTTb) is a reflection of the peak at 8.15 
MeV in Fig. fTlTa) due to this strong correlation and not 
associated with an intermediate state. 

Also possible is the decay ^B* -T-'^He-f-^Lig s. Fig- 
ure [Tljc) shows the reconstructed ^Li excitation-energy 
spectrum determined for the p-a pairs associated with 
the 8.15-MeV peak [gate G7 in Fig. HHa)]). A rough 
background correction was made using gates B3 and BA 
on either side of the 8.15-MeV peak. A broad peak at 
zero excitation energy is observed overlapping strongly 
with the ground state of ^Li. Angular correlations were 
constructed for both the two sequential decay scenar- 
ios, and after background subtractions, are shown in 
Fig. [T31 In this figure 6*12=0 for •^He-f ^Lig decay cor- 
responds to the "^He and a being emitted in the same 
direction. The '^He-^Lig ,5. results is clearly asymmetric 
about cos(0i2)=O and cannot be explained by sequential 
decay. For the p-l-''Be6.73 case, although the distribu- 
tion is symmetric about cos(0i2)=O, we were unable to 
fit it with the range of possible correlations expected for a 
5/2~ intermediate state and thus there is no evidence for 
a sequential decay of this state which suggets the decay 
is 3-body in nature. 
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FIG. 14. (Color online) Background-subtracted angular cor- 
relations associated with possible sequential-decay scenarios 
for the 8.15-MeV state of *B. The 8-15-MeV peak was selected 
with gate G7 in Fig. [11] and the background was estimated 
from B3 U B4. For proton decay to the 6.73-MeV state of 
^Be, the extra gate G5 was applied and for ^He decay to the 
ground state of ^Li, the extra gate G8 was applied. 



2. 2p+ Li exit channel 

The excitation-energy spectrum derived from 2p+^Li 
events in Fig. [15] displays a very prominent peak at 
7.06±0.02 MeV. The excitation energy in Fig. [15] was 
calculated based on the assumption that the ^Li frag- 
ments were created in their ground states and did not 
7 decay. Only one ^Li excited state has a significant 7- 
decay branch, the 3.562-MeV, J''=0+, T=l level whose 
7 branching ratio is practically 100%. Therefore the peak 
observed in Fig.[T5]can correspond to a level of excitation 
energy of either 7.06±0.02 or 10.62±0.02 MeV. 

The intrinsic width of this level is quite narrow as the 
experimental width is similar to the predicted instrumen- 
tal resolution. If we assume a Breit-Wigner line shape, 
we can further constrain the width. The curves in Fig. [15] 
show Monte Carlo predictions including the detector re- 
sponse for intrinsic widths of 50, 75, and 100 keV nor- 
malized to the experimental peak height. The behavior 
in the region of the low-energy tail, where the uncertainty 
due to the background contribution is minimal, suggests 
an intrinsic width of less than 75 keV. 

There are no previously known levels at £'*=7.06 MeV. 
From the mirror nucleus ^Li, we expect only one nar- 
row state near this energy, the mirror of the 6.53-MeV 
,r^A+, r=35 keV level. However, this ^Li state has 
a significant branch to the n+t+a exit channel passing 
through the 4.63-MeV "^Li As the 7.06-MeV state 

is not observed is the mirror exit channel p+^lle+a (see 
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FIG. 15. (Color online) Excitation-energy spectra from 
2p+^Li events. The curves show simulated distributions for 
Breit-Wigner line shapes with the indicated decay widths. 
The lower curve below the data, indicates the background 
which was included in the simulated distributions. 



Fig.llip is seems unlikely that it can be this ,^=4+ state. 

On the other hand, if the peak in Fig. [TSl corresponds 
to a 10.62-MeV state, then we can immediately identify 
it with the known 10.619±0.009 MeV J''=0+ T=2 state 
in ®B with T <60 keV [1^. The almost exact matching of 
the excitation energy and a consistent limit to the width 
strongly support this assignment. 

With this assignment, the ®B state is thus the isobaric 
analog of ^Cg.s. discussed in Sec. IIV Al and its structure 
should be similar. It decays by the emission of two pro- 
tons to the isobaric analog of ^Be^.s. in ^Li. Such a decay 
can conserve isospin only if the two protons are emitted 
promptly in a T=l configuration. Sequential two-proton 
decay does not conserve isospin as there are no energet- 
ically accessible T=3/2 states available in the ''Be inter- 
mediate nucleus (see Fig. llOp . This is an isospin equiv- 
alent to the Goldansky two-proton decay, as the isospin 
allowed intermediate state is energetically forbidden. 

Further elucidation of the nature of this two-proton de- 
cay may be obtained from the correlations between the 
fragments. Unfortunately in this experiment, most ^Li 
fragments were not identified due to a saturation of the 
Si amplifiers, thus most of the 2p+^Li events associated 
with this state were not observed except where the ^Li 
fragment had high kinetic energy (low AE) strongly bi- 
asing any correlation measurement. A future experiment 
is planned to study this correlation. 



3. 2p+d+a exit channel 

The *B isobaric analog state also decays to the 
2p+d+a exit channel. Here, the decay is isospin for- 
bidden and the yield is much smaller than that obtained 
for the 2p+^Li channel (Table |T|, further suggesting that 
the latter is an isospin- allowed decay. The excitation- 
energy spectrum extracted from these events is plotted 



in Fig. [TBla) and shows a narrow peak at the energy of 
the isobaric analog state (arrow). Possible narrow inter- 
mediate states are the ^Be ground state and the 2.18- 
MeV J''=3+ ^Li excited state, although we found no ev- 
idence that the former is associated with this peak and 
in fact used the ^Be ground state as a veto to reduce the 
background. The background was also reduced by gating 
on the reconstructed ^B velocity distribution associated 
with the IAS peak for the 2p-|-^Li events. The spectrum 
labeled "all events" contained these two background re- 
ducing gates. For "all events" , the reconstructed ^Li ex- 
citation spectra from the d-a pairs is shown in Fig. ll6f b) 
where the 2.18-MeV state is quite prominent. Gating on 
this state [gate G4 in Fig.lTUb)]. produces the excitation- 
energy spectrum shown in Fig. I16f a). The peak in this 
case contains only ^50% of the strength of the "all 
events" peak. The origin of the remaining strength is not 
clear, possibly from 3 or 4-body decays, but certainly is 
not associated with sequential decay through a narrow 
state. 

For the strength associated with the 2.18-MeV ^Li 
state, we were not able to identify any ''Be intermediate 
levels. Either the two protons are emitted in an initial 
3-body decay or the ''Be intermediate state(s) are too 
wide to isolate experimentally. Because of the low yields, 
the statistical errors associated with the background- 
subtracted correlations are too large to extract any fur- 
ther information. 



D. 



Be levels 



A high lying *Be state is seen in the p+t+a excitation 
spectrum. The distribution from all detected events is 
displayed in Fig. [TTT a) ( "all events"). A broad peak is 
observed at 22.96 MeV. Such a peak was also observed 
previously following the fragmentation of a ^^Be beam 
and was associated with proton decay to the 4.63-MeV 
J'^=7/2~ ^Li excited state [l2|. This is confirmed in the 
present work. The excitation spectrum obtained from the 
t-a pairs is shown in Fig. [TTT b) where the 4.63-MeV ^Li 
state is quite prominent. Gating on this peak [gate GIO 
inFig.[T7i;b)] produced the "GIO" spectrum in Fig. [TT^a). 
The gate significantly reduced the background, but did 
not reduce the yield in the peak indicating that proton 
decay to the 4.63-MeV state is the dominant decay mode. 

This state has similarities to the 5.93-MeV state in 
®B which proton decays to the 4.57-MeV state in ^Be, 
i.e., the mirror state to the 4.63-MeV in ^Li. Figure [TSl 
compares the angular distributions of 9i2 in the decay of 
these two states. Here the ^B data were scaled to the 
same number of total counts as the ^Be results. Within 
the statistical errors, the correlations are identical sug- 
gesting that the 22.96-MeV state in ®Be is the analog of 
the 5.93-MeV ®B state and thus would have the same pos- 
sible spin assignments (J=3 or 4 ). The intrinsic widths 
of these two states are also comparable; 850±260 and 
680±146 keV for *B and ^Be, respectively. 
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FIG. 16. (Color online) (a) *B excitation-energy spectra ob- 
tained from 2p+d+a events, (b) ®Li excitation-energy spec- 
trum obtained from each d-a pair of each 2p+d+a event. 
The arrows indicate the locations of known states. The gate 
G9 in (b) is used to gate the spectrum in (a). 
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FIG. 17. (Color online) (a) *Be excitation-energy spectra 
from p+t+a events. (b)^Li excitation-energy spectrum from 
the t-a pairs in each p+t+a event. 



In Fig. [Tni we show these levels in an isobar diagram 
of T=l levels. The relative excitation energies of these 
two levels is also consistent with their assignment to the 
same isospin triplet. The corresponding analog in the 
*Li nucleus is not entirely clear. A 6.1-MeV state 
with r ~1 MeV is a possibility; it has been tentatively 
assigned J=3. 



E. 



We observed levels in the p+2a and 2p-j-^Li exit 
channels. Figure BUI shows a level diagram with the levels 
of interest and their decay paths. 



1. p+2a exit channel 

The excitation-energy spectrum derived from p+2a 
events, displayed in Fig. I^ a). shows many peaks. The 
most prominent of these are the ground state and the 
2.345-MeV ^'"=5/2" level. The spectrum above these 
two peaks has been scaled by a factor of three to show 
additional detail. Prominent peaks are also observed at 
11.7 and 14.7 MeV. For the latter, two listed levels can 
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FIG. 18. (Color online) Comparison of angular correlations 
obtained for the 5.92-MeV state in *B and the 22.96-MeV 
state in *Be. The correlations for the *Be state were obtained 
using gates GIO and Gil in Fig. [17] 
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FIG. 19. (Color online) Isobar level diagram of T=\ states 
for A=%. The levels discussed in this work are displayed in 
red. The energy scale on the y axis is with respect to the *Be 
ground state. 
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FIG. 20. (Color online) Level diagram for showing the lev- 
els discussed in this work and their decay paths. The energy 
scale on the y axis is with respect to the ®B ground state. 



contribute, the narrow 14.655-MeV isobaric analog state 
and the wider 14.70-MeV (r=1.35 MeV) state. Judging 
by the shape of the 14.7-MeV peak, it is likely that both 
levels are present. As it is impossible to separate the 
decay modes of these two levels, we will not report on 
any further analysis of this peak. 

The peak at 11.7 MeV is close to a listed 11.640-MeV 
level with a tentative spin assignment of The 
width of this listed level (780±45 keV) is also close to the 
value of 880±80keV extracted from the peak after correc- 
tion for the experimental resolution. For further analysis 
of the correlations associated with this peak, we use the 
gate G12 in Fig. l^lTa) with the regions Bb and BQ on 
either side to estimate the background under the peak. 
The ^Be excitation-energy spectra from the a-a pairs as- 
sociated with each event in the G12 gate is plotted in 
Fig. I^TTb). The "background" was not subtracted from 
this spectrum as the reconstructed *Be [a+a) and ®B 



{p+a+a) excitation energies are highly correlated. This 
is similar to the strong correlation between the '^He-l-a 
and p+'^He-l-ck reconstructed excitation energies found 
in Sec. IIV C II Due to this correlation, the spectra from 
gates Bb and BQ are not good approximations to the 
background associated with gate G12 for the ^Be exci- 
tation energy. In Fig. I^TTb). contributions from the *Be 
ground state and the 3.03-MeV first excited state are 
clearly present. However, the ground-state contribution 
is essentially all background because when we gate on 
this peak and project the ^B excitation-energy distribu- 
tion, no indication of the 11. 7- MeV peak is seen. On the 
other hand, the 3.03 MeV peak is not all background. 
With a similar gating arrangement, we find '^2.8% of 
the strength at 11.7 MeV is be associated with this *Be 
intermediate state. 

Most of the 11.7-MeV ^B peak is thus associated with 
the broad structure above E*{^Be)=b MeV in Fig. CT b). 
This structure overlaps with the wide 11.35-MeV J'^=4+ 
^Be level. The situation here is similar to the 8.15-MeV 
^B state. The ccntroid of the structure in Fig. I^TT b) is 
below the centroid for this level, possibly a consequence 
of sequential feeding. However the *Be and ^B excita- 
tion energies are well correlated and thus the structure 
could just be a consequence of this correlation. The 
alternative intermediate state is ^Li. Fig. I^TTc) shows 
the background-subtracted ^Li excitation-energy spec- 
trum associated with the 11. 7- MeV peak obtained from 
the p~\-a pairs. There are two possible p-\-a pairs from 
each p-\-2a event and both are included in this spectrum. 
A peak corresponding to the ground state of ^Li is clearly 
present indicating that a+^Lig s. is an important decay 
branch. 

The correlations associated with the decay of the 11.7- 
MeV state can be displayed in a manner similar to Fig.fT^ 
where we projected the velocity vectors onto the decay 
plane. To localized the distribution of a particles, the a- 
a relative velocity was made parallel to the Vz axis, and 
as we cannot distinguish the two a particles, we have 
symmetrized the distributions about Vz—0. The veloc- 
ity plot for the 11.7-MeV state is displayed in Fig. I^Wa) 
and can be compared to the corresponding plot for the 
2.345-MeV state in Fig.[22i;c). Unlike Fig. [H] where the 
protons lie on an arc centered at the origin, the protons 
in Fig. I^^ a) lie on two overlapping arcs each centered 
approximately near an a particle. These two arcs indi- 
cate that the protons are not emitted from the ^B center 
of mass, as expected for direct proton decay, but come 
predominantly from the decay of a ^Li intermediate frag- 
ment. Thus the broad peak in the a-a excitation spec- 
trum in Fig. [2lT b) should not be interpreted as an *Be 
excited state. 

The distributions in Figs. l^^ b) and l^^ d) show the re- 
sults of simulations of a sequential disintegration initiated 
by an a+^Lig g. decay. The 2.345-MeV level was previ- 
ously identified as having strong a-f^Lig.s. decay char- 
acter however the two arcs are not separated as the 
two a particles are located too close together. In both 
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simulations we have determined the emission energies in 
the two steps using the R-matrix formahsm as described 
in Sec, mil Both simulations reproduce the experimental 
distribution reasonably well. 

Information about the spin of the 11.7-MeV state can 
be determined from the angular correlations of the two 
decay steps. Unfortunately as the two arcs overlap, we 
cannot always identify which a particle is emitted in the 
first and second steps and thus we cannot uniquely de- 
termine 9i2. Instead we have looked at the angle 92a-p 
defined in Fig. [231 This is the angle between the rela- 
tive velocity vectors of the two a particles {Va~a) and 
the proton velocity (Vp) in the center-of-mass frame. 
Background-subtracted distributions of this quantity are 
shown for both the 11.7 and 2.345-MeV states in Fig. [Ml 
For an initial a-particle decay, unlike proton decay, the 
angular correlations are sensitive to the parity of the ini- 
tial level. Furthermore if we constrain ourselves to a de- 
cays with the lowest possible £ waves, there is no mixing 
to consider and each J'^ value has a unique correlation. 
The prediction obtained for the 2.345-MeV level with its 
known spin of 5/2^ is compared to the experimental data 
in Fig. [MTb). It fits the data extremely well and thus 
again is consistent with sequential decay initiated by an 
a-^Lig.s. decay. In Fig. [Ml^a) for the 11.7-MeV state, we 
show three predictions which fit the experimental data 
with J'^=3/2-, 5/2+, and 7/2". Of these, the 3/2" pre- 
dictions gives the best fit with a x^/t^—^-^- The other two 
cases have x^/z/=2.6 and 3.0 respectively, so J^=3/2^ 
is slightly preferred. We note, that the listed 11.64-McV 
state in ^B has a tentative assignment as J^^7/2^ which 
is the least likely of our three possible fits. 

Figure [Ml shows the energy correlations for the 11.7 
and 2.345-MeV states in ^B. The correlations are pre- 
sented in the Jacobi "Y" system where the "core" in 
Fig. [51 is now the proton and the "protons" in that figure 
are the a particles. The quantity is now the relative 
kinetic energy in the p-a subsystem and, for each event, 
the spectra are incremented twice, once for each the two 
p-a pairs. The curves are the R-matrix predictions us- 
ing parameters for the ^Lig resonance from Ref. [30| . 
These sequential calculations reproduce the data quite 
well, further indicating a sequential character to these 
decays. The strong peak associated with the ^Li ground 
state is clearly resolved in Fig. [SST a) for the 11.7 MeV 
state. 

The discussion of these ^B decays has ignored inter- 
actions between the two decay steps. For the 11.7-MeV 
level, the short life time of the ^Lig.s. intermediate state 
is partially ameliorated by the large kinetic energy re- 
leased in the first step (£^^=10 MeV). With this value 
we obtain Re— 0.12 and ^£=16 fm for this decay. The 
products from the first step thus are separated by '~6 
nuclear radii, on average, before the second steps occurs. 
It is clear that in this case we would expect the sequen- 
tial character of the decay to largely survive final-state 
interactions between the final fragments. 

The sequential character of the 2.345-MeV state is 



more surprising. Here El is much smaller and i?£;=1.88 
and d.E—'i-O fm and this should be a democratic decay. 
The value of d^; is in fact one of the smallest values found 
in this work, even smaller than the value of 6.5 fm ob- 
tained for ^Beg.s. decay. 

Why does this state have such a strong sequential char- 
acter, when ^Bcg.s. decay displays very little? Specifically 
we find that the angular correlations associated with se- 
quential decay are found for the 2.345-MeV ^B state, but 
not for the ^Bcg.s.. A possible reason is that the ^Be^.s. 
decay is approaching a Goldansky decay; most of the 
strength of the intermediate state is energetically inac- 
cessible, whereas all the ^Li^ strength is accessible for 
the 2.345-MeV state. Another consideration is that final 
state interactions between the products would be largest 
for cos(6'2a-p)=±l, i.e., when one of the fragments from 
the second step is directed directly towards the a particle 
from the first step. However, the angular corrections are 
minimal for these cases and so the effect of final state 
interactions will also be minimized. Although the cor- 
relations do maintain a strong sequential character, the 
magnitude of the decay width cannot be explained from 
a R-matrix calculation (Sec. IIIip . The R-matrix predic- 
tion by Barker was a factor of 4 too small f3lj indicating 
the need for a three-body calculation. 

One can also look for ^B levels which decay predomi- 
nantly to ®Be states. The reconstructed ®Be excitation- 
energy spectrum from the a-a pairs of each p+2a event 
is shown in Fig.l^^a). The ^Be ground state and a num- 
ber of excited states indicated by the arrows are clearly 
present. The peak at ^1.5 MeV is not a *Be excited 
state, but rather is associated with the decay of the 2.345- 
MeV ^B excited state. Of particular interest is the strong 
peak associated with the 16.626 and 16.922-MeV mixed 
T=0-|-1 J'^=2~^ states. These are narrow states, but the 
experimental resolution is not sufficient for their separa- 
tion. The peak thus can contain contributions from both 
of these mixed-isospin states. 

The ^B excitation-energy spectrum obtained by gat- 
ing of the E* ~17-MeV peak [gate G13 in Fig. [MJa)] 
is shown in Fig.lMJb). The peak at 16.99±0.04 MeV is 
clearly isolated although it was also present in the un- 
gated spectrum of Fig. [Ufa). This peak could corre- 
spond to a previously known level at 16.71±0.10 MeV or 
the 17.076±0.04-MeV r=3/2 level. Both decays satisfy 
isospin conservation and the former, if it has a signifi- 
cant d-l-^Be branching ratio, may be important in big- 
bang nucleosynthesis [32l | . The magnitude of the yield in 
the observed peak is consistent with that from the un- 
gated spectrum in Fig. l^lT a) indicating that this state 
predominantly decays to one or both of the ^Be T=0-|-1 
levels. Presumably most of the decay is to the 16.626- 
MeV state as only the low-energy tail of the 16.922 MeV 
state is energetically accessible. 

The 16.71-MeV level has only been observed before in 
one experimental study [1^ where it was suggested to be 
the mirror of the 16.67 MeV J''=5/2+ state in ^Be. Ex- 
perimental angular distributions in that work were con- 
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sistent with this spin assignment. The spin of the 17.076- 
MeV state is not hsted in the ENSDF database, however, 
it is probably the mirror of the 16.977- MeV r=3/2 ^Be 
state which has a spin of J^=l/2~. In the two cases we 
would thus expect the first step in the sequential decay 
would be dominated by i—0 and 1 decay, respectively, 
giving rise to different angular correlations. 

We have extracted the angular correlation between 
the first proton emission step and the second 2a decay 
step for the events in the peak located at ~17 MeV in 
Fig. I26r b'). The higher-energy region directly adjacent in 
excitation energy was used for the background subtrac- 
tion. The background-subtracted correlation is plotted 
in Fig. [27] and is compared to calculations for J" =1/2^ 
{£=1) and 7^^=5/2+ (^=0). Of these calculations, the 
data are in much better agreement with the J^=l/2~ 
case, suggesting that the state observed is the T—3/2 
state. The decay to the 16.626 MeV state is Re=0A1 
and (is =45 fm, and the large value of the latter number 
suggests this is approaching a truly sequential decay. 



2. 2p+'^ Li exit channel 

Due to the problem with the saturation of the ampli- 
fiers, only a small fraction of the ''Li fragments could 
be identified. However, there were sufficient statistics 
for the 2p-|-''Li exit channel to isolate a new level in ^B. 
The excitation-energy spectrum determined from 2p4-^Li 
events is displayed in Fig.l^a) and a peak at 20.64±0.10 
MeV is evident with a width of r=450±250 keV. The 
large error is a consequence of an uncertainty of how 
to define the background when fitting this peak. We 
do not know whether the ^Li fragment was created in 
the ground state or the gamma-emitting 0.477-MeV ex- 
cited state. The energy of the level could thus be ei- 
ther 20.64 or 21.12 MeV. We have attempted to see if 
the decay of this level could be associated with a se- 
quential two-proton decay through a *Be excited state. 
As there are two protons, we can reconstruct the two 
possible *Be excitation energies from each 2p-|-^Li event. 
The background-subtracted distribution of both of these 
gated on the 20.64-MeV peak [gate G14 in Fig. [MJa) 
with background gates B7 and B8] is shown in Fig. BST b'). 
Again there is some uncertainty in the background sub- 
traction, but the location of the peak was not affected by 
this. The peak can be identified with the 19.069 MeV, 
r=271 keV, J^=3+ in ^Be. The curves in Fig. ^h) 
show Monte Carlo simulations of the decay through this 
state. The long and short-dashed curves show the con- 
tributions from the correctly and incorrectly identified 
proton associated with *Be decay, and the solid curve 
is the sum of these two. Both the correct and incorrect 
distributions have almost the same average energy, but 
the incorrect distribution is quite wide and the peak po- 
sition is defined by the correctly identified proton which 
has a much narrower distribution. The simulation was 
performed assuming the second decay is to the ground 



state of Li. If we assumed the decay was to 0.477-MeV 
state, then the reconstructed ^Be peak energy would be 
at ~19.5 MeV. No previously known ®Be excited state 
could be found at this energy but it is possible that de- 
cays occurs through an unknown '^Be excited state, how- 
ever it seems more likely that the decay is to the ^Li 
ground state through the known 19.07-MeV state. 



F. "'C levels 

Information on ^'^C states produced by neutron pickup 
and more complicated reactions was obtained from the 
2p+2a exit channel. This is the only channel available 
for particle-decaying states with E* < 15 MeV. The 
excitation-energy spectrum from all of these events is dis- 
played in Fig. BW a) and is dominated by a peak at 9.69 
MeV, although other smaller peaks are evident. In our 
previous work, levels for this exit channel were revealed 
by gating on the presence of various narrow intermedi- 
ate states [H]. Figures l^b). [^c). and [HJd) show 
the excitation-energy distribution gated on the ^B 
ground state, the ^Be ground state, and the 2.345-MeV 
^B excited state, respectively. 

The arrows in Fig. I^Wb) show the locations of two 
states (5.222 and 6.553 MeV) that were excited via the 
inelastic scattering of a beam and which decay by 
proton emission to the ^B ground state j23i] . These states 
dominated the total excitation energy spectrum in 
inelastic scattering, but here they are weakly populated 
and require a gate by the ^^g.s. to enhance them [see 
Fig.[29lb)]. As these states are expected to have 2a clus- 
ter structure (23j , then it is not surprising that a neutron- 
pickup by (which does not have any 2a cluster struc- 
ture) does not populate these states strongly. 

The widths of these states also appear to be wider 
(by 50-70%) than we would expect based on the intrinsic 
widths given in Ref. f23| and the simulated experimental 
resolution. Possibly these peaks are actually doublets 
which would not be surprising as other levels near these 
energies are expected based on the mirror nucleus ^Be. 
At higher excitation energies in Fig. I29f b). there may 
also be indications of other states, though the statistical 
errors are large. We note the absence of any significant 
contribution from the 9.69-MeV level in this spectrum. 

A 6.56-MeV state which decayed via the a-l-^Be^.s. 
channel was also observed in the inelastic excitation of a 
beam [2^, [s^l- Its location is indicated by the arrow 
in Fig. BW c). This state is also expected to have strong 
cluster structure and is clearly only weakly populated in 
the present experiment. In Fig. [29lc). the spectrum is 
again dominated by the 9.69-MeV peak, indicating that 
a-t-^Bcg.s. is an important decay branch for this level. 

The spectrum in Fig. I^^ d) is also dominated by the 
9.69-MeV level indicating that this level has a second de- 
cay branch (p-l-^B2.345). Also a smaller peak appears at 
8.5 MeV. The arrow in this figure shows the location of 
an 8.4 MeV state that was previously shown to decay 
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through the 2.345-MeV intermediate level. However, 
the intrinsic width of that state was ~1 MeV, signifi- 
cantly wider than the present peak. The experimental 
width of the peak in this work ('--^300 keV) is consistent 
with the simulated resolution. Within the statistical un- 
certainties we estimate F <200 keV. 

Other small peaks are also suggested in the total dis- 
tribution in Fig. ^a). The peaks at 10.48 and 11.44 
MeV, indicated by the arrows in this figure, are previ- 
ously unknown. Their widths are consistent with the 
experimental resolution. The decay of these states is not 
known; they do not feature prominently in the distribu- 
tions gated by the long-live intermediates. 



1. 9.69-MeV state 

The strong 9.69-MeV state is wider than the other ob- 
served states; from the simulation we estimate its 
intrinsic width is 490 keV. It was not observed via ^"^C 
inelastic scattering in Refs. [H, H^l suggesting that it 
does not have a strong cluster structure and maybe is 
more single-particle in nature. In any case we will show 
its decay modes are quite interesting. In the following 
we give an analysis of correlations associated with this 
peak using gate CIS and background gates 59 and BIO 
in Fig. I^^ a) . All the distributions discussed in the re- 
mainding of this section are background subtracted. 

As we have observed in Fig. [MT c) and [Mtd), the 
9.69 MeV level has strength for decays into the a-l-^Beg 
and p-l-^B2.345 channels. Estimates of the branching ra- 
tios for these two channels are based on the ^Be and ^B 
excitation-energy spectra gated on the 9.69-MeV state 
displayed in Fig. [301 One can clearly see peaks associ- 
ated with ^Be ground state [Fig. Wi^)] and the 2.345- 
MeV state [Fig. in ^B. The yield associated with 

the ground state of ^B is clearly minimal. The yield in 
the ^Be and ^B2.345 peaks were fit with Gaussians with 
smooth backgrounds. From these yields and correcting 
for the detector efhciency, we determine branching ra- 
tios of 35% and 17% for a-f ^Bcg.s. and p-l-^B2.345 decay 
paths, respectively. Both of these decay branches are 
safely sequential with (i_B=136 and 200 fm, respectively. 
The determination of the decay mechanism of the re- 
maining ~50% of the yield is more difficult and, to isolate 
this component, we have vetoed events associated with 
the ^Bcg s. and ^B2.435 peaks in the following analysis. 

The distribution of Epp, the relative kinetic energy 
between the two protons, is plotted in Fig. [3TJ This 
distribution shows a strong enhancement at small rel- 
ative energies [Epp <1.0 MeV), the diproton region, but 
also displays a tail to significantly higher values. Fig- 
ures [32r a) and IHSl b) show the distributions of ^Be and 
^Li excitation energy determined from the detected a- 
a pairs and all 4 combinations of p-a pairs. The total 
distributions (circular data points) have strong overlaps 
with the wide 3.03-MeV ^Be and the ^Li ground-state 
resonances, respectively. The distributions are rather in- 



sensitive to whether the events are associated with the 
diproton peak in Fig. [3T] or associated with the higher- 
energy tail. To illustrate this, the squared data points 
in Fig. [5Wa) and [5Wb) show the distributions gated on 
the diproton peak {Epp <1.0 MeV). Within the statisti- 
cal uncertainties, they are the same shape as the ungated 
distributions. 

The strong overlap with the ^lAg ^. resonance may 
initially suggest that the ^"^C fragment decays via the 
^Lig ,,..-|-^Lig channel. However in such a case, the ex- 
citation energy reconstructed from an a particle and a 
proton originating from different ^Li fragments should 
populate the ^Li excitation energies significantly higher 
than ^Lig s. resonance region. However, the experimental 
data indicated that all four a-p pairs are consistent with 
the ^Lig.s. resonance. The decay thus appears to be four- 
body in nature and cannot be described by a sequential 
process. A possible interpretation is that all four p-a 
pairs are simultaneously in ^Lig ,5. resonances and some 
fraction of the p-p pairs are tightly correlated in a " dipro- 
ton". Possibly, the a-a pair may also be in a 3.03-MeV 
^Be resonance. On the other hand, the overlap with the 
^Lig s. and other resonances may just be coincidental. If 
all of these associations with resonances are correct, it is 
not clear that a configuration exists which would satisfy 
all of the implied angular momentum couplings. To dis- 
play the correlations graphically, we note that each p+2a 
triplet in each 2p~\'2a event defines a plane in velocity 
space. With the two protons, we define two such planes 
in each event. The distribution of A(/)pp, the angle be- 
tween these two planes is plotted in Fig. [33la). This dis- 
tribution resembles the relative-proton-energy distribu- 
tion in Fig. [3T] with a diproton region and a high-energy 
tail, and in fact the quantities Epp and A0pp are strongly 
correlated. 

In order to help visualize the correlations between each 
proton and the two alpha particles, we have employed 
the procedure used to generate Fig. [T^ For each p+2a 
triplet, the a-particle and proton velocity vectors are pro- 
jected into their decay plane. Subsequently they are ro- 
tated in this plane to align the locations of the a parti- 
cles. The 2-dimensional velocity distributions are plotted 
in Fig. [33] for gates on the diproton region {A(j)pp < 40°) 
and the high-energy tail {A(j)pp > 40°). 

Both distributions are quite similar and thus again 
there is little difference in the p-a correlations for the 
diproton and high-energy tail regions of the p-p relative 
energy. This is further illustrated in Fig. I33f b) where 
the projected Vx distributions for these two regions are 
compared. Again we see very little difference in the re- 
sults for the two gates. The average location of the two 
a particles and the proton form an isosceles triangle in 
their defining plane; i.e., the two p-a average relative ve- 
locities are identical and thus all p-a pairs can have the 
same reconstructed ^Li excitation energy. 
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V. DISCUSSION AND CONCLUSIONS 

This work highlights the wide range of many-particle 
decays of ground and excited states in light nuclei. It 
is not clear that there is a simple way of the predicting 
when the decays should be sequential or prompt. Cer- 
tainly the presence of possible narrow intermediate states 
does not guarantee that the system will take advantage 
of them. For example the 9.69-MeV state in has 
^Beg.s. and ^62.34 intermediate states with T ~90 keV; 
yet 50% of the time it bypasses these and undergoes a 
4-body decay. The 6.56 MeV state also has com- 
peting prompt and sequential decay branches (23j . The 
nature of the decay is certainly related to the structure of 
the state and, in principle with suitable models, this in- 
formation can be used to improve structure calculations. 
Presently there are only a few cases where prompt 3- 
body decay has been modeled in a sophisticated manner, 
as for example the quantum-mechanical cluster calcula- 
tion where comparisons to experimental correlations 
are useful. The prompt 4-body decay poses an even more 
difficult theoretical problem. 

In this work we have studied many-particle exit chan- 
nels associated with states formed following the interac- 
tions of an E/A^70 MeV beam with a ^Be target. 
The decay products were detected in the HiRA array 
and the states were identified from the invariant mass of 
the detected fragments. Correlations between these frag- 
ments were used to deduce the nature of the decay and 
search for long-lived intermediate states associated with 
sequential decay processes. For such decays, angular cor- 
relations between the fragments provide information on 
the spin of the level. 

The ground state of *C disintegrates into four pro- 
tons and an a-particle. This decay was found to proceed 
through the ^Be^.^. intermediate state, with two steps of 
2-proton decay. The correlations between the protons in 
the first step exhibit a significant enhancement in a region 
of phase space where the two protons have small relative 



energy. The correlations in the second step were found to 
be consistent with ^Be^.s. 2-proton decay obtained from 
an independent measurement using a ''Be beam. 

Evidence suggests that the isobaric analog of this state 
in ^B, undergoes 2-proton decay to the isobaric analog of 
^Bcg s. in ^Li. This is the first case of two-proton decay 
between two isobaric analog states. The exact nature 
of the 2-proton decay awaits future measurements, but 
it may constitute a Goldansky-type decay with a single 
proton emission being energetically allowed but isospin 
forbidden. 

A 9.69-MeV state of ^°C was found, which -50% of the 
time, undergoes a prompt 4-body decay into the 2p+2a 
channel. There is a strong diproton-like correlation be- 
tween the two protons, at the same time, all four pos- 
sible p-a pairs have energies consistent with a ^Lig 
resonance. No sequential decay scenario is possible to 
produce 4 simultaneous ^Lig s. resonances. 

Evidence for a prompt 3-body decay of the 8.15-MeV 
state of ^B was found. On the other hand, a number of 
states with three-body exit channels were observed that 
are well described as sequential decays through interme- 
diate states. These include the isobaric pair, ^Bs.ga and 
8Be22.9, and the 11.7, 16.9 and 20.6 MeV states of ^B. 
The correlations in the decay of the 2.345-MeV state of 
^B to the p + 2a channel were found to have strong 
sequential character associated with a ^Li intermediate, 
even though the lifetime of this intermediate state is very 
short and, at the very least, we would expect there to be 
very important final-state interactions. 

Finally, new mass excesses and decay widths were ob- 
tained for the ground states of ^B and *C. 
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FIG. 21. (Color online) (a) Excitation-energy spectra ol 
®B derived from p+2a events. (b) The *Be excitation- 
energy spectra obtained from the a-a pairs for events in 
the 11.7-MeV (G12 gate), (c) Background subtracted ^Li 
excitation-energy spectra for events in the 11.7-MeV ^B peak. 
The arrows indicate the locations of known levels. 



^' (a)E*=il.7MeV, ex'p 




„ 2 

cn 

c 

E 
o 



(b) E*=il.7 MeV, Sim 

'iLC^JJA 


i 







(c) E*=: 

a 


1.345 rlfleV, e'xp 




■ 



(d) E*= 


^.345 lifleV, e'xp 




1 



-2 2 4 
[cm/ns] 



-2 2 4 
[cm/ns] 



FIG. 22. (Color online) Distributions of the velocities of the 
three fragments for the p+2a exit channel associated with 
(a,b) the 11.7-MeV and (c,d) the 2.345-MeV states of ^B. The 
distributions are projected on the plane of the decay with the 
two a-particle velocities localized. The experimental distri- 
butions in (a,c) are compared to simulated distributions for 
a-^Lig.s. decay in (b,d). 




FIG. 23. (Color online) Schematic showing the definition of 
the angle 92a-p- 
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FIG. 24. (Color online) Background-subtracted angular cor- 
rections measured in the decay of the (a) 11.7-MeV and (b) 
the 2.345-MeV states of ^B. The curves indicate calculated 
distributions assuming an initial a-^Lig,s. decay. In (a), the 
calculated results are given for a number of assumed spins for 
the 11.7-MeV state. 
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FIG. 25. (Color online) Background-subtracted energy cor- 
rections in the Jacobi "Y" system measured in the decay of the 
(a) 11.7-MeV and (b) the 2.345-MeV states of ''B. The curves 
indicate calculated distributions assuming a-^Lig.s. decay in 
the R-matrix approximation. 
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FIG. 26. (Color online) (a) Excitation-energy distribution for 
*Be fragments reconstructed from the a-a pairs of all p+2a 
events, (b) The excitation-energy distributions obtained 
with the G13 gate on the *Be excitation energy. 
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FIG. 27. (Color online) Angular correlation for the 16.99 MeV 
state in ®B. Data points indicate the background-subtracted 
experimental results and the curves shows the results ex- 
pected for the indicated level spins. 
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FIG. 28. (Color online) (a) Excitation-energy distribution for 
®B fragments reconstructed from 2p-|-^Li events, (b) The *Be 
excitation-energy distributions obtained with the Gl gate on 
the ^Be excitation energy. The distributions contains contri- 
butions from both possible p-^Li combinations. 
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FIG. 29. (Color online) Excitation-energy distributions for 
^"C fragments reconstructed from 2p+2a events. Panel (a) is 
for all events, while (b), (c) and (d) are gated on the presence 
of ^Bg.s., '^Beg.s., and ^Ba. 345 intermediate states. Arrows 
show the locations of the peaks discussed in this work. The 
gate G15 used to select the 9.69-MeV peak and the gates B9 
and -BIO used for background subtraction are indicted. 
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FIG. 30. (Color online) Background-subtracted excitation- 
energy distributions for (a) ''B and (b) ®Be fragments pro- 
duced in the decay of the 9.69-MeV level of ^"C determined 
from p+2a and 2p+a triplets. 
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FIG. 31. (Color online) Distribution of relative kinetic energy 
between the two protons for the 9.69-MeV level in ^'^C. 
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FIG. 32. (Color online) Background-subtracted excitation- 
energy distributions for (a) *Be and (b) ^Li fragments pro- 
duced in the decay of the 9.69-MeV level of ^"C determined 
from a-a and p-a pairs. For the latter case, the distribu- 
tion contains all four possible p-a pairs associated with each 
detected 2p-\-2a event. The distributions indicated by the 
circular data points have no restriction on the p-p relative 
energy, while for the squared-shaped data points Epp < 1.0 
MeV. The arrows with error bars show the mean energy and 
width of the indicated resonances. 
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FIG. 33. (Color online) Background-subtracted distributions 
associated with the decay of the 9.69-MeV state in ^°C. (a) 
Distribution of the angle A^pp. This is the angle between the 
two planes in velocity space defined by the two a particles 
and each of the protons, (b) For the two indicated gates 
on A(^pp, the distributions of the x component of the proton 
velocity (see text). For comparison purposes, the results for 
A<jf>pp > 40° have been normalized to the same number of 
counts as in the A0pp < 40° distribution. 
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FIG. 34. (Color online) Velocity distributions of a particles 
and protons produced in the decay of the 9.69-MeV state in 
^"C. For each 2a+p triplet in each detected 2Q+2p event, 
the velocities have been projected onto the plane of the three 
fragments and rotated in this plane to align the locations of 
the a particles. The results are shown for (a) A(/)pp < 40° 
and (b) A(j!)pp > 40° where the A(/)pp is the angle between the 
two planes defined by the two protons. 



